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Abstract 

By the extension of the 3-dimensional analytical solutions of Goldreich and Weber [6] with 
adiabatic exponent 7 = 4/3, to the (classical) Euler-Poisson equations without cosmological 
constant, the self-similar (almost re-collapsing) time-periodic solutions with negative cosmo- 
logical constant (A < 0) are constructed. The solutions with time-periodicity are novel. On 
basing these solutions, the time-periodic and almost re-collapsing model is conjectured, for 
some gaseous stars. 
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1 Introduction 

The evolution of a self-gravitating fluid, such as gaseous stars, can be formulated by the isentropic 
Euler-Poisson systems of the following form: 

( p t +V • (pu) =0, 

(pu) t +V • (pit® u) + VP =-/>V$, (1) 
A$(t,x) = a(N)p — A, 

where a(N) is a constant related to the unit ball in R N : a (2) — 2ir and a(3) = An. And as 
usual, p = p{t, x) and u = uit, x) € R w are the density, the velocity respectively. P = Pip) is the 
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pressure function. While A is the cosmological constant, for A > 0, the space is open; for A < 0, 
the space is closed; for A = the space is flat. 

In the above systems, the self-gravitational potential field $ = <&(<, x) is determined by the 
density p itself, through the Poisson equation (jTJ)3 . For A 7 ^ = 3, (JT|) are the classical (non-relativistic) 
descriptions of a galaxy, in astrophysics. See [T], 0, [7], for details about the systems. 

The 7-law can be applied on the pressure term P(p), i.e. 

P(p) = Kp-r = ^, (2) 

7 

which is a common hypothesis. The constant 7 = cp/c v > 1, where cp, c v are the specific heats 
per unit mass under constant pressure and constant volume respectively, is the ratio of the specific 
heats, that is, the adiabatic exponent in @. In particular, the fluid is called isothermal if 7 = 1. 
Historically in astrophysics for zero constant (A = 0) background, for the 3-dimensional space, the 
hydrostatic equilibrium specified by u = 0, was studied. According to [2], the ratio between the 
core density p(0) and the mean density 'p for 6/5 < 7 < 2 is given by 

= (—y(z)) , (3) 



p(o) V » 



z=z 



where y is the solution of the Lane-Emden equation with n = 1/(7 — 1) for 7 > 1 : 

y(z) + -y(z) + y(z) n = 0, y(0) = a > 0, y(0) = 0, n = (4) 
z 7 — 1 

and zq is the first zero of y(zo) = 0. We can solve the Lane-Emden equation analytically only for 

l--z 2 , n = 0; 
.6 
sin z 

n = 1; (5) 



Uanal (^) 



1, n = 5, 



I V/1 + ZV3 

and for the other values, only numerical values can be obtained. It can be shown that for n < 5, 
the radius of polytropic models is finite; for n > 5, the radius is infinite. 

As we are interested in the solutions in radial symmetry, the Poisson equation ^ 3 is transformed 
to 

r N - l <S> rr (t, x) + (N- 1) r N - 2 $ T =u (N) (p - A)^" 1 , (6) 
*r = ?3© [ r (p(t,s)-A)s N - 1 ds. (7) 



For the analytical collapsing solutions, the well-known results are constructed by Goldreich and 
Weber [6] in 1980 for 7 = 4/3 in R 3 . The analytical collapsing solutions without compact support 
for 7 = 1 in R 2 are constructed by Yuen pj] in 2008. 

For simplicity, we take the constant A = — for A < 0. For the background knowledge for the 
Euler-Poisson equations with cosmological constant A, the interested readers may refer [5], [9]. 
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Our main results are able to obtain the analytically time-periodic solutions for the Euler-Poisson 
equations with the negative cosmological constant (A < 0) in spherical symmetry: 



u = -V(t,r), (8) 
r 

1/2 

where the diameter r :— f 5Zf=i x i) m -^ 3 > 

2 

Pt + Vp r + P V r + - P V = 0, 

r (9) 
P (V t + VV r ) + -§^(Kp i /' i ) = --^j- fo(p(t, s) + -^)s 2 ds, 

with the extension of the analytical solutions of Goldreich and Weber [5] in the theorem: 

Theorem 1 For the 3- dimensional Euler-Poisson equations in spherical symmetry (0|), there exists 
a family of solutions, for 7 — 4/3: 



,(*,r)=(* y( =fr )8 ' f0rr<am - ,V(t,r)=^r; 



0, fora(t)Z v <r. 



a W = Z(tF ~ a (°) = a o > °' a (°)= ai ! 
y(z)+|y(z) + f y(z) 3 ^, y(0) = a > 0, y(0) = 0, 



(10) 



where K > 0, p = 3A/(4if), and the finite Z v is the first zero point ofy(z). 
In particular, 

(1) A = and a± = 0, the solutions collapse in the finite time t = ir/2; 

(2) A < 0, the solutions are non-trivially time-periodic, except the case with ao — \/A and a\ = 0.; 

(3) A > and a\ < 0, the solutions collapse in a finite time T . 

2 Time-Periodic Solutions 

In this section, before presenting the proof of Theorem [T] we prepare some lemmas. First, we 
obtain a general class of solutions for the mass equation in spherical symmetry ©i- 

Lemma 2 For the 3- dimensional equation of conservation of mass in spherical symmetry: 

Pt + V Pr + pV r + - P V = 0, (11) 

r 



there exist solutions, 



«-'>-W- VM -$r- (12 > 



with the form f > e C 1 and a(t) > £ C 1 . 
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Proof. We just plug (O into fTTjl. Then 



Pt + V Pr + P V r + -pV (13) 



r 



-3a(t)f(r/a(t)) a(t)r f (r / a(t)) 



(14) 



(15) 



a{ty a{tf 

a(t)rf(r/a(t j) f (r / a(t)) a(t) 2 f (r / a(t)) a(t) 
a{t) a(i) 4 a{t) 3 a(t) r a(tf a(t) 

= 0. (16) 

The proof is completed. ■ 

The following core lemma is needed to show the cyclic phenomena of the solutions (|T0|) . 

Lemma 3 For the ordinary differential equation, 

a(t) - — X — + a{t) = 0, a(0) = a > 0, a(0) = a u (17) 
a(t) 

with A > 0, there exist non-trivially periodic solutions except the case with clq = s/A and a\ = 0. 

Proof. To the equation (fl7| . multiply a(t) and then integrate it: 

«W 2 , A a{tf _ 



with the constant 9=^- + ^ + ^>0. 
We define the kinetic energy: 



F km := (19) 



and the potential energy: 



The total energy is conserved: 



j t (F ki n+F pot )=0. (21) 

Observe the potential energy has only one global minimum at a m ; n = S^A, for a(t) £ (0,+oo). 
Therefore, by the classical energy method for conservative systems (in Section 4.3 of [8]), the 
solutions have the closed trajectory. We calculate the time for traveling the closed orbit: 

t = f dt= 2 r ~ da{t) - 2 Mt) (22) 

Jo y w J 2[9 _ (_*_ + A ml „ ^ 2[fl _ ( _A J + £ffl!) ] 

where a m j n = fnf (a(t)) and a max = sup(et(t)). 

*>o t>0 

Let G(f) := G := 0-( z -A T _+ i-^) > and G x := e-f^+i*^) > 0. 
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Except the case with ao = \f\ and a\ = 0, the equation (f2"2")l becomes 



T 



a m i n +e 



2da(t) 



2da(t) 



a(t) 2 



2da(t) 



2^-(^T + ^)] 



< 



sup 

a min^ a (£)^ a min + ^ 



sup 

( — e<a(f)<a n 



sup 



A a(t) 



Go 



X «(*) 



a(t)- 
1 



V2dG(t) 
V2dG{t) 



2da(t) 



a m i n +e 



w - a(t) 



v Co 

"7T 



2da(t) 



sup 

-e<a(£) <a n 



A «(*) 



V2 



< oo . 



(23) 
(24) 

(25) 

(26) 
(27) 

(28) 
(29) 



Therefore, the solutions to the differential equation (fiTf are time-periodic. 
The proof is completed ■ 

In the solutions ([T0| . due to the properties of the modified Emden equation (flQ|) 9. the collapsing 
phenomena are observed. The following lemma can be proved immediately. 



Lemma 4 For the ordinary differential equation, 

X 



a(t) 



a{tf 



a{f), < o(0) = a , a(0) = oi < 0, 



(30) 



with A > 0, 



there exists a finite time T such that lim aft) = 0. 

Proof. If the claim is not true, for all t > 0, we have a(t) > 0. 
But the ordinary differential equation (|30|) becomes 

■1 



a (ty- 



a{tf 



-a(t) < 0. 



(31) 



(32) 



(33) 



After a sufficient large time, there exists a finite time T such that lim a(t) — 0. A contraction is 

t->T- 

met. 

The proof is completed. ■ 



This shows d(t) is decreasing function: d(t) < d(ti) for all t > t-y > 0, 

a{t) < o(ti) < 0. 

Thus, the solution is bounded by 

a(t) = / d(s)ds + a(ti) < d(ti)t + ao- 
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Here we are ready to give proof of Theorem [TJ 
Proof of Theorem [TJ From Lemma [51 we can easily check that (|10p satisfy ©i- We plug the 
solutions (fTO)) . into the momentum equation 



p{V t + VV r ) + K%-pV* + ^ / (p(t, s) - A)s 2 ds 



a(t) 



Or 



1 + —f / (p(t, s) - A)s 2 ds 



a{tf 



a{ty 



+ 



a(t) 
Xr AttA 



+ AK viM viwj) top 



3/ s n 2 , 47rA P 
27 ( ^) )S ^ ~ 



P 



a(t)3 



-pr + 4if p 



a(i) 3 a(t) 2 r 2 a(t) 2 
u 

o(t)2 + r 2 a{tfj V ( a(t) } 3 P 



a(*) 2 



r 



r/a(t) 



y(s) 6 s 2 ds 



a{t) 2 ^ 

Here, we use the property of a(t): 



and denote 



... . —A AttA — A 

a(<) = ^ + — a(<) = ^- a(<) ' 



7* 4"7r 
Q(— ) := Q(.t) = -Xx + 4Ky(x)+— / y(s)V<fa. 
a(t) x z ' 



Differentiate Q(x) with respect to x, 



(34) 
(35) 

s 2 ds (36) 
(37) 
(38) 

(39) 
(40) 

(41) 
(42) 



Q(x) = -A + 4Ky(x)+4iry(x) 3 



2-47T 



y(s) 6 s^ds 



X 

Xx + 4Ky(x) - Kax + / y(s) 3 s 2 ds 
x 2 J 



-Q(x). 



(43) 

(44) 
(45) 



With lim Q(x) = Q(0) = 0, this implies that Q(x) = 0. 

The above result is due to the fact that we choose the Lane-Emden equation: 

y( Z ) + ly( Z ) + §y(z) 3 =p., M = 

y(0) = a > 0, y(0)= 0. 



(46) 
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We have shown that there exists the family of the solutions with compact support, in Theorem [TJ 
with the well-known results about the Lane-Emden equation (|4"u| [2] [7]. 

On the other hand, for A = — §tt < 0, it is clear for A = in the theorem is true. For the linear 
ordinary differential equation: 

a(t) + a(t) = 0, a(0) = a > 0, d(0) = 0, (47) 

the solution is 

a(t) = ao cost. (48) 

The function a(?) achieves zero and the solutions collapse in the finite time tt/2. 

By using Lemma [3] about a(t), we have , for A < 0, the periodic solutions are non-trivial except 

the case with ao = s/A and a% = 0. 

By using Lemma |4j we have, for A > and a\ < the solutions collapse in a finite time T. 
This completes the proof. ■ 

Remark 5 If we consider the system with frictional damping, 

p t +V • [pu) =0, 

(pu) t +V • (pu <g> u) + /3pu + VP =-pV3>, (49) 
A$(t,f) = a(A r ) j o-A, 

wii/i t/ie constant j3 > 0, 

i/ie corresponding ordinary differential equation for the solutions are: 

a(t) + pa(t)=-^r - a(t), a(0) = a > 0, a(0)= oi, (50) 

m 

For A < 0, the function a(t) oscillates around the minimum a m i n = v^A of the potential energy 
except the case with ao = v^A and a\ = 0. And it is asymptotically stable. 

Remark 6 In general, for the Euler-Piosson equations in R , it is clear to have the corresponding 
solutions: 

for N = 2 and 7=1, 

p(t,r) = 4 m z v{r,a[t) \V(t,r)=^r- 

o(t)=^-r - o(t), a(0) = a > 0, a(0)= a i; (51) 
a(t) 

1 27T V ^ 

y(z)+-y(z)+—e = p, j/(0) = a, y(0) = 0, 
z K 

where p = 2X/K, fTJ^; 



8 



M.W.Yuen 



for N > 3 and 7 = (2N - 2)/N, 

1 



Ar)= l^)" ,M -'«"<«, p(t , r) 4, 

0, for a(t)Z„ < r. 



a(t)-- 

2/(2)+ 



2/(z) = M> 2/(0) = a > 0, 2/(0) = 0, 



(52) 



(2JV-2)Jf 

w/iere fi = [N(N - 2)A]/(2AT - 2)K and the finite is the first zero ofy(z) J$. 

With X < 0, the solutions\51\ and\5^ are time-periodic for the negative cosmological constant. 



Remark 7 Under some initial conditions for A > 0, the collapsing solutions are obtained. The 
case is similar to the system without cosmological constant. 



3 Discussion 

By the separation method, the solutions are constructed in its close relatives (the Euler-Poisson 
equations with frictional damping, Navier-Stokes, Euler equations and the ones with frictional 
damping term) [TU] and [Hj. However, we emphasis on that it is novel to have the time-periodic 
pattern in the solution structure with A < 0. 

On the other hand, in the statement (2) of the theorem, with some initial conditions, for 
example, < a(0) = e << 1 and d(0) = 0, the density becomes 

a 3 

p(0,0) = — »M 0) (53) 

where Mo is an arbitrary constant. 

It may provide the possible evolutionary model with A < constant in mathematical theory, that 
the universe expands and then almost re-collapses (The density at the origin can be greater than 
any given constant p(T, 0) >> M ), in the finite time periodically. The solutions may be like a 
black hole (in the sense, p(T,0) — +00 with a finite time T). The black hole seems to be the end 
of the evolution of the periodic solutions. However actually, there exists no authentic black hole 
in our particular time-periodic solutions. 

As the time-periodic effect is due to the negative cosmological constant in our model, more 
solutions with this pattern are expected for other 7 values. Then further works may be done by 
the computing simulations about the stability of the numerical solutions. If some gaseous stars (a 
galaxy) obey the 7- law (7 = 4/3), it may provide an alternative explanation about that, the time- 
periodic solutions coincide with the expansion segment (the red-shift effect) in a short time. We 
notice that it is the significant difference between the traditional re-collapsing model with negative 
cosmological constant A < [3]. 
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